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We report on measurements performed on an apparatus aimed to study the Casimir force in the 
cylinder-plane configuration. The electrostatic calibrations evidence anomalous behaviors in the 
dependence of the electrostatic force and the minimizing potential upon distance. We discuss analo- 
gies and differences of these anomalies with respect to those already observed in the sphere-plane 
configuration. At the smallest explored distances we observe frequency shifts of non-Coulombian 
nature preventing the measurement of the Casimir force in the same range. We also report on mea- 
surements performed in the parallel plane configuration, showing that the dependence on distance 
of the minimizing potential, if present at all, is milder than in the sphere-plane or cylinder-plane 
geometries. General considerations on the interplay between the distance-dependent minimizing po- 
C**) , tential and the precision of Casimir force measurements in the range relevant to detect the thermal 

corrections for all geometries are finally reported. 
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PACS numbers: 12.20.Fv, 03.70.+k, 04.80.Cc, ll.10.Wx 



I. INTRODUCTION 



Casimir forces [l| have been investigated since their inception as a macroscopic test of the irreducible fluctuations 
associated to quantum fields. They are geometrical in character, as they originate from the possibility to confine 
and shape the energy density of quantum fluctuations using proper boundary conditions. So far, the experimental 
£NJ ■ attention has been mainly focused on the original parallel plate configuration 0, E| , and the sphere- plane geometry 
[H-Q, apart from the only experiment performed in a crossed-cylinder configuration Further motivations to 

. pursue precision studies of the Casimir forces are related to the possibility to discover new forces of strength similar 
or larger than the Newtonian gravitational coupling but with short-range or different distance scaling and expected to 
act, according to various models, in the submillimctcr range fill ]. This can be considered part of a broader program 
aimed at testing deviations from Newtonian gravitation in the nonrelativistic limit (l2l . [l~3| . In the micrometer range 
the dominant source of background to non- Newtonian gravitational forces is provided by Casimir forces [lil [ToT ] (see 
also [H, [TtJ for reviews) and to discover such forces, or at least to provide reliable limits on their existence, one must 
control at the highest level of accuracy all systematic sources of deviation from the idealized case analyzed by Casimir 
in his original paper. Major sources of systematic errors that are still considered under partial control are the use of 
the proximity force approximation (PFA) [3 [l9| for curved surfaces, the presence of electric forces not reflected in 
the purely Coulombian contribution, such as patch effects (2fjj - [2^ , and the combined effect of finite conductivity and 
' finite temperature (25l - [27| . 

The use of the PFA has been discussed at length in the literature, with several alternative methods developed to 
overcome its limitations, and exact solutions have been found in particular curved geometries. It is generally assumed 
that the PFA for the Casimir force differs from the exact result by an amount smaller than 0.1 %, an assumption 
compatible with the results of a dedicated experiment [HI . The control on the PFA used to assess limits to Yukawian 
non-Newtonian gravity has not been addressed until very recently, although it has been used for many years [p, [29l - 
HH . It has been argued that the usual form of PFA cannot be extended unambiguously to volumetric forces [35| . and 



thereafter an alternative form of the PFA has been discussed [36| , which, however, has been shown in [37[ to coincide 
with the exact formula for geometries in which one of the two bodies has translational invariancc. 

The presence of electric forces not incorporated in the Coulombian contribution has been discussed extensively in 
the literature, in particular in [2(| . Anomalies in the electrostatic calibrations of the sphere-plane configuration have 
been evidenced for large radii of curvature of the sphere and small distances from the planar surface [38l . |39| . This 
has triggered discussions about the nature and the universality of the observed anomalies [H, H(J , and the situation 
is still far from being clarified. The anomalous exponent optimizing the fit of the electrostatic calibrations in [3^, [39| 
has not been found in another experiment using spheres of much smaller diameter located at similar distances from 
the planar surface [4l| , while the dependence of the minimizing potential on the sphere-plane separation reported in 
[H| has been confirmed in [4l| , and for crystalline Ge in [42j (see also [39| for a discussion of some unpublished data 
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from former experiments). The spatial and temporal variabilities of the minimizing potential has been evidenced in a 
centimeter-size torsional balance [43j |. which confirms the necessity for a detailed knowledge of the surfaces and their 
preparation p3 - l47j . 

The finite temperature contribution added to the quantum fluctuations has originated a lengthy debate about the 
interplay of the thermal contribution with the finite conductivity properties of the surfaces (see for instance [48l — 
l59l | for the initial steps of the debate). On the experimental side, attempts to evidence the thermal contribution 
discriminating various models have been reported for the sphere-plane geometry while proposals using torsional 
balances in the parallel plane configuration [60l l6lj are under development. A dedicated experiment in the parallel 
plane configuration using microresonators |62| and a low-frequency heterodyne technique [63} has been limited so far 
from patch charges 64 1 . 

Recently, the use of a configuration with intermediate features between the parallel plate and the sphere-plane 
ones, i.e. the cylinder-plane geometry, has been proposed, and its experimental feasibility was investigated at gaps 
of the order of 20 Lim, limited by the roughness of the metallic surfaces (65l - [67| . This geometry is very relevant 
from the theoretical viewpoint, since an exact solution for the Casimir force has been found IrM 6 9 f . also providing 
another example of curved geometry in which PFA may be tested against numerical techniques [70L l7l| . In this paper 
we report on the results of electrostatic calibrations for an apparatus using the cylinder-plane geometry in a range 
of distances of relevance for measuring the Casimir force. We observe a background originating frequency shifts of 
amplitude large enough to overwhelm the downshift expected from the Casimir force. We also discuss the distance 
dependence of the contact potential in both the cylinder-plane and parallel plates configurations. The minimizing 
potential shows no significant distance dependence in the parallel plates configuration with respect to corresponding 
cases of the cylinder-plane and sphere-plane configurations. 

The paper is organized as follows: in Sec, II we briefly recall the Coulomb force in the cylinder-plane geometry, and 
report on the upgrades to the apparatus with respect to the one described in (66| , its overall sensitivity performance, 
its geometrical characterization, and the parallclization technique. In Sec. Ill we describe results from electrostatic 
calibrations, showing that in analogy to the sphere-plane case we observe that (a) the optimal exponent for fitting 
the dependence of the Coulomb coupling on distance is not the one expected from the idealized situation, at least at 
the smallest explored distances between the cylinder and the plane and (b) that the minimizing potential depends on 
distance. We then describe the data analysis leading to force residuals after subtraction of the Coulombian contri- 
bution. At the smallest explored gaps we observe residual frequency shifts of amplitude large enough to prevent the 
measurement of the Casimir force. The presence of shifts neither of Coulomb nor of Casimir origin has been confirmed 
by implementing a measurement strategy consisting in progressively approaching the two surfaces at constant bias 
voltage, and measuring the resonator frequency. In Sec. IV we discuss possible explanations for the unexpected 




FIG. 1: (Color online) (Left) Schematic view of the experimental set-up with two lateral piezoelectric actuators (PZTl and 
PZTr) for fine adjustment of parallelization, and a central piezoelectric actuator (PZTc) for controlling the cylinder-plane 
distance. On the top of the central actuator a platform is located with screws to adjust the horizontal position of the Au coated 
cylindrical lens (C). The resonator (R) is on the top of the cylindrical lens, and below the optical fiber (F) which is attached 
to a piezoelectric actuator (PZTf). (Right) Close-up image of the cylindrical lens and resonator region. The cylindrical lens 
has a radius of curvature of a — 12 mm and a length L — 4 mm, smaller than the width of the resonator of 10 mm. 
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FIG. 2: (Color online) Assessment of the parallelization through capacitance measurements. The cylinder-plane capacitance is 
shown vs. the difference between the steps traveled by the left and right actuators in a differential mode. Each data point is 
the average of 20 measurements with an integration time of 1 s. The error bar represents the mean standard deviation on each 
data set. 

scaling law of the Coulomb interaction with distance. We then describe in Sec. V electrostatic calibrations taken in 
a parallel plane geometry aimed at evidencing the dependence on distance of the minimizing potential also in this 
configuration, thereby completing the analysis for the three most common geometries of experimental interest. The 
relevance of measuring and modeling the dependence of the minimizing potential on distance in order to detect the 
thermal contribution to the Casimir force is discussed in Sec. VI. In the conclusion we put our findings in the more 
general framework of the recent observations of systematic effects highlighting possible future developments for the 
cylinder-plane geometry 

II. CYLINDER-PLANE CONFIGURATION: GEOMETRICAL CONSIDERATIONS 

The theory related to the cylinder-plane geometry and a description of the apparatus have been the subject of a 
former paper [r36j , while further details of the measurement technique have been reported in (38l . l39l . l67l . [72j . The 
main difference from the former tests is due to the use of high-quality Au-coated cylindrical lenses, with planarity 
and roughness comparable to the ones of the resonator, therefore allowing us to reach submicrometer gaps. 

The electrostatic force between a conducting cylinder (of length L, radius a. with L 3> a) parallel to a conducting 
planar surface, separated by a gap d and kept at a fixed electrostatic potential difference V is [z|[ 



where A = yW — aP and h = d + a, with the approximate expression in the right-hand side (r.h.s.) valid in the limit 
d«a, also coinciding with the result expected from the PFA for electrostatics. Since both the cantilever and the 
cylinder have finite size, the length L in the previous expressions should be replaced by an effective length L e g that 
characterizes the relative exposure between the cantilever and the lens, (i.e. the minimum between the width of the 
cantilever and the length of the cylinder). In our initial experimental attempts this corresponded to the cantilever 
width. However, by visual and optical microscope analysis we noticed the presence of sharp irregularities at the 
border of the cantilever, most likely originated by the laser cutting process of the Si wafer. We have then chosen 
cylindrical lenses of length L = 4 mm smaller than the resonator width of 10 mm (see Fig. 1), as the lenses seem to 
have more regular borders, as visible at the optical microscope. At this point we should note that the curvature of our 
cylindrical lens is a = 12 mm, larger than its length L. This implies that border effects in the electrostatic interaction 
between the cylinder and the cantilever can be important, and the exact logarithmic expression of the electrostatic 
force in Eq.((T]) should not hold in our configuration. However, as long as the conditions for the PFA hold (d <C a) 
the approximate expression for the force as given by the r.h.s. in Eq.([T]) should apply, irrespective of the relative 
magnitude of L and a. Other important geometrical issues are the possible nonperfect parallelization between the 
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cylinder and the cantilever and the use of a cylindrical lens rather than a full cylinder. In [66[, the correction to the 
PFA expression for the force [r.h.s. of Eq. (JTJ] due to nonparallelism was computed 
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where a = L sm9/2d, 9 is the deviation angle from ideal parallelism, and in this nonparallel case the distance d 
between the cantilever and the cylinder is measured from the midpoint along the axis of the cylinder. The use of 
a cylindrical lens rather than a full cylinder can be simply evaluated in PFA resulting in a subleading PFA 
correction (of the order of 5 x 10 -2 for a typical separation of d/a = 10~ 3 ), and will be then discarded in what follows. 

To calibrate the apparatus, a controllable electrostatic force is generated by applying bias voltages between the 
cantilever and the cylindrical lens. At a given separation, the frequency of a resonant mode characterized by an 
effective mass m e g is measured both with (v) and without [vq) the presence of the voltage. This allows for the 



evaluation of the square frequency difference Av 



related to the voltage V in the PFA (d <IC a) as 
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where Vq is the minimizing potential. For small tilting angles a <C 1 the correction to the frequency shift due 
to nonparallelism has a quadratic dependence on a, given as 1 + 35a 2 /24. The parallelization procedure is one- 
dimensional, thereby simpler than in the case of two flat surfaces. Rather than measuring the frequency shift induced 
by a constant bias voltage as discussed in [66j . we have opted to monitor the cylinder-plane capacitance using a 
capacitivc AC bridge [72j]. Using the PFA, the capacitance between the cylinder and the plane is given by 



_ 27reo£effVQ (V 1 + a - \/l - a) 
V2d a 

The tilting angle 9 (and therefore a in turn) is controlled through two motorized actuators acting in differential 
mode, with the goal to keep the average distance d (the separation between the plate and the cylinder as measured 
from the midpoint of the cylinder) constant while changing the tilting angle (see Fig. 1). The precision of the 
achieved parallelism is then determined by the quality of the fitting of the capacitance versus differential steps 
number, propagated to determine the dispersion on the number of steps at minima compatible with the fitting error. 

The precision can be improved by minimizing stray capacitance between various contacts since the precision of 
the capacitance meter is some percentage of the total capacitance in the system, usually 0.05%. An example of 
parallelization through capacitance measurements is shown in Fig. [2 The use of long-range actuators has been 
shown to be problematic due to large hysteresis, and therefore we have opted for a fine tuning with the use of two 
piezoelectric actuators as in the schematics shown in the left plot of Fig. I. An optical microscope allows to monitor 
the quality and cleanness of the two surfaces and for a visual, qualitative assessment of the parallelization. An a 
posteriori, off-line check is obtained by fitting the electrostatic curves of the frequency shift using Eq. which 
results in a value of 9 compatible with zero within a fitting error of 59 = 10~ 3 radians. 



III. ELECTROSTATIC CALIBRATIONS AND RESIDUALS ANALYSIS IN THE CYLINDER-PLANE 

GEOMETRY 

An important prerequisite to any Casimir force measurement is the execution of high-quality electrostatic calibra- 
tions. Forces are always indirectly measured via their functional relationship to more accessible observables, such as 
detection of deflection angles, voltages required to keep the apparatus at rest in closed loop schemes, or shifts of the 
frequency of a mode of a resonator as in our case. It is then crucial to convert such quantities more directly observed 
into the corresponding force signal by means of well-known and controllable physical signals, for instance by applying 
external bias voltages and comparing the measured forces with the Coulombian interactions between macroscopic 
conducting bodies. 

The electrostatic calibration starts by finding the best parallelization condition at a given nominal cylinder-plane 
separation d using the capacitance technique described previously. Then the parallelization is further fine tuned by 
adjusting the two lateral piezoelectric actuators PZTl.r. Once the optimal parallel condition is obtained, PZTl,r are 
left untouched, and the separation d is changed via the central piezoelectric actuator PZTq. This procedure keeps 
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the parallclization at the optimal value. As seen previously, in the perfect parallel situation the frequency shift of the 
cantilever due to pure electrostatics takes the form 



2 _ 3e 0x /^L cff (V - V ) 2 _ 2 

^ "WW - - Kcliv Vo) ■ (5) 

The displacement of the central piezo PZTc depends linearly on the voltage applied VpzT, and therefore the absolute 
gap is given by d = f3(Vp ZT — Vpzt), where /? = (91.9 ± 0.9) nm/V is the actuation coefficient of the piezoelectric 
transducer, and Vp ZT is the PZT voltage required to make contact between the two surfaces. At a given distance, the 
electrostatic calibration has been performed by measuring the frequency shift induced by a range of electric voltages 
V applied between the two surfaces. The curvature coefficient K e \ = 3eoy/aL e g/16\/2Trm e f{d 5 ^ 2 and the minimizing 
potential Vq can then be obtained by fitting the data with 



"d = "*i(V, V PZT ) = 4 ~ K ci (V PZT ) x (V - V ) 2 , (6) 

^el(VpZT) = 7(^PZT - V PZT )- 5/2 , (7) 

where 7 = 3eo^/aL e ff/16\/27rm c f[P 5 / 2 . This fitting procedure allows the determination of the absolute distance d 
once the fitting parameter Vp ZT is obtained, and the measurement of the contact potential Vq as a function of d. A 
typical data plot of K e \ vs. Vpzt is shown in the left plot of Fig. [31 The blue curves are the best fits using Eq. ©, 
and they deviate significantly from the data points, both including or excluding a curvature offset representing 
a hypothetical background electric field. Moreover, the effective mass calculated from the fitting parameter is 30- 
50 times larger than the physical mass, much larger than the expected value of the effective mass which should 
be comparable or smaller than the physical mass of the resonator. If instead the power exponent is left as a free 
parameter, rather than being fixed at 2.5, a new fitting curve with the exponent in the 0.9 to 1.3 range is obtained 
(the red curves in Fig. [3]) . This deviation from the expected exponent for the Coulombian force has been confirmed 
to exist in all our electrostatic calibrations data. 

Table I shows the fitting parameters of the electrostatic calibrations for five runs both when the exponent is fixed 
and left as a free parameter. The effective mass for the latter case is not well defined because of the deviation of the 
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FIG. 3: (Color online) (Left) Plots of curvature coefficient K e \ vs. Vpzt in cylindrical-plane electrostatic measurements obtained 
with the curvature technique, and best fits with the expected Coulombian interaction having a 2.5 exponent (blue/dark gray 
continuous curve), and with a power-law functional dependence in which the exponent is instead a free parameter (red/light 
gray continuous curve). The dashed curves are obtained by also constraining to zero an offset for K c \, representing a possible 
curvature present even at a large distance between the cylinder and the plane, for instance due to stray environmental electric 
fields. (Right) Plot of the reduced x 2 obtained by dividing x 2 by the degrees of freedom (DOF) vs. the value of the power 
exponent in the case of the offset = 0. Data are taken as described in [3Sl . [3^ ] , with a modification of the acquisition code 
for faster data acquisition and smaller uncertainty from drift effects. In the routine used for data acquisition in [38I. |39|]. the 
bias voltage always returns to V after each measurement. In the new routine instead, it changes from +14 to —14 in steps of 
SV where 14 and 5V (i.e. the maximum bias voltage and its minimum step of variation during the calibration) are specified in 
advance in the code. This allows to double the acquisition speed for a targeted span of voltage values. Although this sacrifices 
some accuracies in the subsequent fitting of the data, a shorter data acquisition time is preferable considering the relatively 
large drift experienced during the entire duration of a typical run. 
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Run 1 


Run 2 


Run 3 


Run 4 


Run 5 


fixed exponent 
V P ° ZT (V) 
d min (nm) 

mes (g) 

■J 2 /nm? 
X /L)(ji> 


2.5 
79.52±0.06 
590 
6.63 
fin 


2.5 
73.08±0.06 
504 
9.25 

1 1 fi 
110 


2.5 
56.76±0.03 
491 
8.87 

Zoo 


2.5 
71.05±0.02 
464 
9.46 

Olo 


2.5 
65.84±0.02 
477 
9.01 

R7fi 
o ( D 


free exponent 
Vp° ZT (V) 
dmin (nm) 

x 2 /dof 


1.30±0.03 
74.37±0.04 

116 

5.3 


0.97±0.01 
68.07±0.02 

43 

11 


0.93±0.01 
51.78±0.01 

33 

19 


1.00±0.01 
66.45±0.01 

42 

42 


0.94±0.01 
61.03±0.01 

34 

31 



TABLE I: Fitting parameters of the electrostatic calibrations using the curvature technique. 



exponent from 2.5, thus it is not listed. As shown in Table I, the exponent, when left as a free parameter, is always 
smaller than the theoretical value of 2.5, with a reduced % 2 smaller by about an order of magnitude with respect to 
the one expected from the Coulombian scaling. The relatively large value of x 2 also indicates that the errors may 
be underestimated, although this does not affect our conclusions about the relative comparison between Coulombian 
and optimal exponents. While small deviations from 2.5 are expected considering all the less than ideal conditions 
such as imperfect parallclization and thermal and mechanical drifts, such a significant difference (an average value of 
1.03 versus 2.5) cannot be explained as small deviations from ideality. 

Electrostatic calibrations have also been performed with an alternative technique consisting in directly measuring 
the resonance frequency as the separation gap is decreased maintaining a constant bias voltage. This so-called fast- 
approach measurement technique has the advantage of a faster data acquisition, resulting in a mitigation of the long- 
term thermal or mechanical drifts, and provides an alternative to check the distance dependence of the electrostatic 
force. As shown in the left plot of Fig. 21 the optimal exponent obtained from the minimum of the reduced \ 2 curve 
is around 0.89. Although this technique seems to be slightly less sensitive to the exponent, as shown from the softer 
dependence of \ 2 (the value of x 2 for the exponent of 2.5 being only three times larger than the minimum x 2 obtained 
at the exponent of 0.89), the optimal value of the exponent is consistent with the results from the previous technique 
based upon electrostatic calibrations. The fact that a significant deviation of the exponent from 2.5 is still observed 
even in fast-approach measurements limits or rules out the possibility that systematic effects, such as artifacts from 
fitting the parabolic dependence in Eq. (|5|), or long-term drifts such as thermal expansions or relaxation of the PZT 
actuator, may be responsible for this anomalous behavior. 

As can be seen in Table I, the smallest distance achieved is of the order of 500 nm based on the fitting with a 
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FIG. 4: (Color online) (Left) Squared resonance frequency versus Vpzt for the fast-approach calibration technique. The data 
were obtained using a measurement protocol in which a constant bias voltage equal to 4 V was applied across the two surfaces 
and the resonance frequency was progressively measured from the farthest to the closest distances. The dashed blue curve is the 
fitting with a 2.5 exponent while the continuous red curve is obtained leaving the exponent as a free parameter, whose optimal 
value turned out to be about 0.89. In the inset the reduced x 2 is plotted versus the value of the free exponent. (Right) Test of 
the electrostatic scaling law for the cylinder-plane geometry at large distances with the curvature technique. Plot of curvature 
coefficient K c \ versus Vpzt. The red curve is the fit with 2.5 exponent, the blue is the fit with the optimal exponent 1.84. In 
the inset the reduced x 2 is plotted versus the value of the free exponent, with a minimum of x 2 obtained for an exponent equal 
to 1.84. 
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FIG. 5: (Color online) Optimal exponents of five electrostatic calibrations versus the distance of the closest data point used in 
the fitting, with the continuous line representing a spline curve obtained by considering all data points, weighted by their error 
bars. 



fixed exponent of 2.5. In a previous measurement using wider cylindrical lenses with a larger radius of curvature 
the smallest distance reached was around 1 /im, and no significant deviation of power exponent from 2.5 has been 
observed. This suggests that the exceptionally small exponent may be a result of the smaller gaps reachable with the 
cylindrical lenses of smaller width and radius of curvature. Although the absolute distances obtained from the fitting 
with exponent 2.5 cannot be fully trusted in light of the relatively inaccurate fitting, they can be still considered as 
reliable enough to estimate the gap separation. 

Electrostatic calibrations with the same lens were performed at relatively large distances, as shown in the right 
plot of Fig. 4. In the data presented in Table I, the explored distances ranged from about 500 nm to 7.3 /im. The 
new calibrations were instead performed in the range of about 5.1 to 21.6 fjxci. This obviously results in much smaller 
frequency shifts and larger error bars in the values of K e \, yet it is evident from the fitting that an exponent of 2.5 
is more satisfactory for these large-distance data. Furthermore, the fitting parameters with an exponent fixed at 2.5 
gives an effective mass of 1.2 g which is much closer to the estimated physical mass of the resonator of 0.2 g. If the 
power exponent is instead left as a free parameter, a value of 1.84±0.06 is obtained. The large-distance data span 
smaller ranges of K e \ and therefore are less sensitive to the power exponent and as shown in the inset on the left plot 
of Fig. HI the difference between the fitting curve of 2.5 exponent and that of 1.84 exponent is not very significant 
(corresponding to a reduced x 2 of 2.3 versus 4.9). Nevertheless, the fitting makes evident that the optimal exponent 
at large distance is almost two times larger than the one obtained from the small-distance data. We also investigated 
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FIG. 6: (Color online) Dependence of the minimizing potential on distance in the cylinder-plane geometry. (Left) Minimizing 
potential Vq versus Vpzt from a typical electrostatic calibration measurement. (Right) Minimizing potential Vo versus distance 
d from various runs. 
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FIG. 7: Plots of the curvature coefficient K e \ (left) and the minimizing potential Vo (right) versus time at constant Vpzt- The 
inset of the right plot is the histogram of Vo which shows that Vo roughly follows a Gaussian distribution centered around 0.163 
V with a full width at half maximum of 0.023 V (standard deviation 0.01 V). 

how the optimal exponent changes if a subset of data, rather than the entire set, is used in the fitting procedure. 
In particular, we have fitted subsets of data obtained by progressively removing points at the smallest distances. As 
shown in Fig. [5l all five runs share the same trend showing that the optimal exponent increases when the number of 
the removed point of closest distance used in the fitting increases. The absolute distances in the plot were obtained 
from the fitting with fixed 2.5 exponent. One feature which is noticeable in the figure is a relatively sharp increase 
in the value of the optimal exponent for a distance range of 500 to 600 nm. This could explain why such a large 
deviation of the optimal exponent from 2.5 has not been observed in our earlier measurements with a larger radius 
of curvature cylinders in which we managed to reach minimum distances of only about 1 /zm. We have also noticed 
a strong correlation between the value of the effective mass m c fj and the exponent. The effective mass obtained from 
the fitting increases when the distance between the two surfaces decreases. The dependence on distance of the optimal 
exponent is in contrast to the case of the sphere-plane measurements in which a relatively constant optimal exponent 
was observed uniformly over the entire range of explored distance [HI, [3i| . All fits are performed with a weight equal 
to l/crf , where <7j is the standard deviation of K e \ from the parabola fitting. 

In our experiment we have also studied the possible distance variability of the minimizing potential Vo {i.e. the 
voltage difference between the cylindrical and the planar conducting surfaces which is minimizing the electrostatic 
force, as described for instance in [38l . |39(). At larger gaps we have observed an approximate linear relationship 
between the residual potential Vo and Vpzt in the cylinder-plane configuration. With closer approaching (larger 
values of Vpzt), Vo tends to have a nearly flat dependence on the distance, as is visible in Fig. [51 both on a single 
run (left plot) and on various runs obtained in different days (right plot). 

To estimate the effect of long-term drifts we have studied the time dependence of the electrostatic curvature 
coefficient K c \ and the minimizing potential Vo, without nominally changing the cylinder-plane gap distance, as 
shown in Fig. [7] The curvature coefficient shows temporal variations of order 50 %, while the minimizing potential 
Vo does not show any evident dependency, rather it fluctuates in a relatively small range of values. 

From the electrostatic calibration measurements it is also possible to extract the electrostatic-free frequency v as 
can be seen from Eq. ©. In principle, when the distance between the two surfaces is small enough, one should expect 
a downward shift of this residual frequency due to attractive Casimir forces. However, as can be seen in Fig. |8l no 
significant downshift of vo has been observed, rather a sharp upshift was instead observed in one case. This upshift 
indicates a very strong repulsive force between the two surfaces at small distances, which could be from either a 
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184 
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368 


Vbias 


= 3V 


-0.91±0.02 


-1.52±0.10 


-1.56±0.17 


-1.91±0.32 


-2.05±0.51 


Vbias 


= 4V 


-0.89±0.03 


-1.20±0.08 


-1.73±0.19 


-1.89±0.30 


-1.83±0.36 


Vbias 


= 5V 


-0.98±0.02 


-1.70±0.04 


-1.60±0.06 


-1.61±0.10 


-1.62±0.16 



TABLE II: Optimal exponents for various values of minimum distance used in the fitting of the data from fast approach 
measurements. 
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FIG. 8: (Color online) Residuals from electrostatic calibrations with parabola method. (Left) Plot of residual frequencies from 
various runs of electrostatic calibration measurements vs. distance. Because of the dependence of the resonant frequency on 
temperature, different runs have different frequencies even at the largest gaps at which no residual force is expected. (Right) 
Same plot but with a common baseline chosen in such a way that the data points at the farthest distance have a common 
central value, with v s — uq — fo(Gf max ). 



repulsive component of the Van der Waals forces or simply a soft contact of the surfaces. To better understand the 
short-distance behavior and to reduce the effect of thermal drifts, we have performed a series of measurements using 
the fast-approach technique mentioned previously. 

As shown in Fig.|9l with a constant bias voltage Vb; as =4 V, the frequency of the resonator decreases as the cylindrical 
lens approaches the resonator, as qualitatively expected for the attractive Coulomb force, then followed by a sharp 
increase, most likely from contact of the surfaces. If the Coulomb force is the dominant force, then the downward 
part of the curve can be fitted using Eq. [5j As mentioned earlier, a fixed exponent of 2.5 produces a marginal fit while 
the optimal exponent is around 0.9. However, we should keep in mind that the exponent is not supposed to be 2.5 
in the first place if at small distances some other forces (either the Casimir force, patch forces, or corrections to the 
standard Coulomb interaction mentioned in the previous section) become large enough to compete with the Coulomb 
force. Based on this remark, we further analyze the data by removing data points at the smallest distances, to test 
the stability of the fit with the Coulomb force. In these conditions the optimal exponent increases, as shown in Table 
II. In analogy to the previous analysis on electrostatic calibration data, although the optimal exponent is still smaller 
than 2.5, the value is approaching 2.5 within the relatively large error bars. This indicates once again that the sources 
of discrepancy between the data and the expected Coulomb force are localized at the smallest distances. 
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FIG. 9: (Left) Plot of frequency versus Vpzt with a constant bias voltage Vbias = 4 V obtained using the fast-approach method. 
(Right) Plot of frequency versus Vpzt with a constant bias voltage Vbias = 0.15 V approximately compensating the value of 
the minimizing potential expected at the smaller gaps. No evidence for downshifts attributable to charge-independent forces is 
visible until the frequency increases due to direct contact between the two surfaces. 
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FIG. 10: (Color online) Residuals to the frequency shifts with the fast-approach technique. (Left) Plots of residual frequency 
square (y 1 - u^ t ) for a constant bias voltage Vbi as = 4 V versus Vpzt and different values of the maximum value of Vpzt used for 
the data analysis. (Right) Plot of the corresponding force versus distance. Both the absolute force and the absolute distance 
are inferred by the Coulomb fitting at larger distances. 



We then can investigate another kind of residual of the data to analyze the effect of data taken at the smallest 
distances. Instead of relaxing the exponent finding its optimal value, the exponent is kept fixed at 2.5 when fitting 
the larger distance portion of the data, and the residuals at small distances are evaluated. This residual analysis is 
shown in Fig. [10] for one run, with the soft contact occurring around Vpzt = 88.8 V. If the Coulomb force is the only 
dominant force for the whole downshifted part before the soft contact, then the fitting curve obtained when excluding 
a small region of data prior to contact, for instance all the data corresponding to a PZT voltage larger than Vpzt = 
87.8 V, should be able to predict the data obtained by excluding distances corresponding to a further Volt of Vpzt 
removed, with residuals centered around zero. However, by doing so there is clearly a nonzero downshift residual 
before the soft contact takes over. In principle this residual could come from fitting artifacts, in particular it could 
depend on the interval chosen for fitting the data with the Coulomb force. However, when more points were removed, 
this nonzero residual appears to be stabilized. Considering that the frequency shifts with this fast-approach technique 
may capture all possible forces acting on the resonator, this electrostatic residual analysis indicates that there are 
forces other than the expected contribution from the applied constant bias voltage that caused a further downshift 
in the frequency. To test whether the residual force is correlated to the external bias electric field, measurements 
were also performed with Vbias = 0.15 V, corresponding to the average value of the residual potential Vq at small 
distances as shown in Fig. [6j No noticeable downshift was observed, as seen in the right plot of Fig. [9] Furthermore, 




FIG. 11: (Color online) Fast-approach measurements with various values of the bias voltage. Plot of the resonator frequency 
versus distance with bias voltage Vbias of 3 V (black circle points), 4 V (red square points), and 5 V (green diamond points). 
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FIG. 12: (Color online) Plots of residual frequency square (y - i/ fit ) from the fast-approach measurement with Vbi as = 3 V 
(black circle points), 4 V (red square points) and 5 V (green diamond points) versus Vpzt (left) and distance (right). The 
fittings were done after removing the closest 184 nm data. 



to rule out possible changes of configurations (such as parallelism or distance drifts) in runs taken in different days, 
fast-approach measurements were performed with different Vbi as applied within the same run, (i.e., at each position 
three values of T4>ias were applied and the respective frequencies measured, as shown in Fig. [TT] with bias voltages of 
3, 4, and 5 V). This confirms that larger bias voltages result in larger frequency shifts at the same distance. The 
residual analyses were carried out and the results with the fitting after removing the closest 184 nm data are shown 
in Fig. 1121 The extra downshifts were present in all three curves, with peak value approximately quadratic in the 
external bias voltage. 

Although this residual force F rcs is attractive, it cannot be identified with the sought Casimir force. Apart from the 
absence of a comparable signal in the residual frequency analysis with the electrostatic calibration technique (Fig. [8]), 
Fros seems to be dependent on the applied bias voltage as shown in Fig. [12] Moreover, best fits of this residual force 
with power-law expressions indicate that it is required an exponent which is much larger than that of the Casimir 
force. The fact that this extra force depends on the applied bias voltage indicates that it may also be present in 
electrostatic calibrations and could lead to an anomalous exponent. 

It should be noticed that within our statistics of runs performed with the fast-approach technique (about 20) we 
have also observed a couple of runs in which residuals gave rise to short-distance upshifts. This could be explained by 
the presence of anomalous distance drifts due to external factors such as the environmental temperature. To assess 
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FIG. 13: (Color online) Plots of the resonator frequency versus time in the case of large separation (black upper points), and 
small separation (red lower points), the latter obtained by imposing an external bias voltage of 4 V at a nominally constant 
distance. While the first plot reflects the intrinsic resonator change, for instance due to temperature drifts or internal creeps, 
the second also includes the effect of drifts over time in the cylinder-plane separation distance. 
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FIG. 14: Optimal exponent ds. the distance of the closest data point used in the fitting of the capacitance vs. distance curve for 
the numerical simulation of the finite size cylindrical lens in front of the finite size planar resonator, the geometry corresponding 
to the actual experimental setup as shown in Fig. 1. 



this effect and how much it affects the amplitude of the observed frequency shifts in the residuals, we show in Fig. [13] 
the typical amplitude of the frequency fluctuations. This allows for disentangling the intrinsic drifts due to changes 
in the resonator frequency, obtained by monitoring the resonator with a distance from the cylinder large enough to 
make negligible its influence, and drifts due to changes in the cylinder-plane distance, with a bias voltage of 4 V 
intermediate between the two extreme values of voltages applied in the fast-approach measurements. 



IV. POSSIBLE EXPLANATIONS FOR THE ANOMALOUS EXPONENT 



In the following section we discuss possible causes of the anomalous exponent obtained in our electrostatic calibra- 
tions of the cylinder-plane geometry. In particular, we consider edge effects, local deformations from the idealized 
geometry, electric forces with steeper distance scaling than the Coulomb force, and electrostatic patch effects. 



A. Edge effects 



One possible reason for a deviation from the ideal Coulomb prediction of 2.5 for the exponent of the frequency 
shift versus distance is the fact that a finite size cylindrical lens was used in the experiment, instead of a long whole 
cylinder. In this situation edge effects may not be negligible, and Eq. ([T]) would not be a good approximation to 
the actual electrostatic frequency shift. While edge effects have been discussed for Casimir forces with the general 
world- line approach in [74| . we have not found former discussions of this effect in the electrostatic calibrations in 
Casimir experiments. Using the COMSOL numerical package, we have conducted numerical simulations in which the 
precise geometry of the measurements was used, and the capacitance between the cylindrical lens and the resonator 
was evaluated at different distances [lf|. By neglecting edge effects, the power exponent for the capacitance versus 
distance would be 0.5. We repeated the same analysis as shown in Fig. [5] with these data, and Fig. [F4l shows how 
the optimal exponent changes if data points at the smallest distances are progressively removed from the fitting. The 
optimal exponent obtained for the capacitance from the COMSOL simulation for the finite-size geometry deviates 
more from the ideal value of 0.5 as the cylinder-plane separation is increased (i.e. the regime in which edge effects 
are more pronounced). In contrast, the optimal exponent obtained experimentally in our electrostatic calibrations 
deviates more strongly from the ideal value as the separation is decreased. Therefore, border effects cannot explain, 
either quantitatively or qualitatively, the anomalous exponent observed in our measurements. 



B. Local geometrical deformations 



Another possible reason for the deviation of the exponent with respect to the ideal value 2.5 is the presence of 
geometrical deformations in the shape of the cylinder. As the experiment is performed with a very large cylinder 
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FIG. 15: Two examples of deformations of the cylindrical surfaces, 
on the right, the cylinder has a tip (not to scale). 



On the left we show a cylindrical lens with a flat deformation, 



(a=12mm), the surface may present local deformations at the submillimeter scale, which could induce strong deviations 
in the exponent. A similar point has been raised for the case of a large sphere in front of a plane in [40(, where it 
was remarked that departures from the ideal spherical surface may noticeably affect the exponent of the electrostatic 
calibration. 

To illustrate this point, we have computed the electrostatic force between a deformed cylinder and a plane using the 
PFA. Let us first assume that, in the region of minimum distance between surfaces, the cylinder has a flat deformation 
of width 2b (see left plot in Fig. [T5|) . For simplicity we assume that the same deformation is all along the length L of 
the cylinder. In this case it is possible to show that 



e L e g(V~V ) 2 d 2 



4tt 2 



m c fi 



dd 2 



[fmc(d- 



2a 



) + f PP (d)} 



where 



Jinc (a) = y —j- arctan W — 



(8) 



(9) 



is the contribution of the (incomplete) cylinder, and f pp (d) = b/d is the contribution of the flat deformation. Let us 
assume that Ai^j = —A/d B . Depending on the relation between the size b of the deformation and the range of distances 
d of the calibration, we expect that the exponent B will interpolate between the ideal value 2.5 (cylinder-plane) at 
relatively large distances, and 3 (parallel plates) at short distances (although the interpolation is not necessarily a 
monotonic function). When b = 10 2 /xm, and the fit is performed for d between 0.5 and 2 fim, the exponent becomes 
B = 2.8, bigger than that of the ideal cylinder-plane geometry. Then this kind of deformation does not help to explain 
the observed anomalous exponent. 

On the other hand, if the cylinder has a deformation with the form of a tip (see right plot in Fig. I15[) . the exponents 
are, in general, considerably smaller than the ideal one, and can explain at least part of the anomaly. Indeed, let us 
assume that the deformation consists of a triangular tip of width 2b and height b' . In this case, PFA gives 



e L eS (V-V ) 2 d 2 



4tt 2 



m cf [ 



where 



/1 



tip 



dd 2 



In 1 



[finc(d + b') + f tip (d)}, 



(10) 



(11) 



If the height of the tip is much larger than d, the contribution of the incomplete cylinder is almost irrelevant, since 
the cylinder is shifted upward and its electrostatic energy becomes almost independent of d in this regime. The 
main contribution comes from the tip, and the mild logarithmic dependence of the energy with the distance implies 
an exponent of around B « 2. This can be easily confirmed by performing fits of Eq. (|10j) . For example, for the 
same parameters as previously, with b' = 6, we obtain an exponent B = 2.0. Sharper tips may produce even smaller 
exponents, although the PFA becomes unreliable for very thin tips. 

The main conclusion of the PFA estimations is that, as expected, deformations of the cylindrical surface may change 
appreciably the exponent 2.5 of the electrostatic calibration. Although the previous examples do not explain the full 
discrepancy between the ideal prediction and the experimental data, this is certainly a crucial point to be taken into 
account in future experiments. 
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C. Additional electric forces 



In the residuals analysis of electrostatic calibrations described in the previous section, a residual force which seems 
to depend on the applied bias voltage was observed. The presence of such a force could lead to an anomalous exponent. 
Let us consider a hypothetical scenario in which the square of frequency shift A^ s for a cylinder-plane geometry has 
the following dependence on distance: 



a 2 ( a \ 



(12) 



in which p > 2.5. This means that besides the expected electrostatic interaction between the two surfaces, there is 
another electric force which follows a higher power law upon the cylinder-plane separation. Let us generate a set of 
pseudo-data following Eq. (fT2"j). and let us try to fit the curvature coefficient K p = aid -2 5 + a 2 d~ p with: 



A" 



a(d — do) 



(13) 



where a, do and q are fitting parameters. In Figs. [16] and [Ijj we show the results of the fitting using different 
values of as/ax versus the distance of the closest point used in the fitting with p = 5. Not surprisingly, the optimal 
value of q from the fitting approaches p at small distances, and approaches 2.5 at large distances. However, it is 
not monotonically decreasing from p to 2.5 when the distance increases, but first goes below 2.5 and then slowly 
approaches 2.5. The 0:2/0:1 = 5 and 10 curves in Fig. [TBI resemble the results from our electrostatic calibrations data 
as shown in Fig. Therefore it is possible that an extra force with a power law dependence on distance steeper than 
the Coulombian one is responsible for the strong deviation of the optimal exponent from 2.5. The value of a from 
the fitting also suffers from similar issues. In the electrostatic calibration measurement, the exponent is fixed at 2.5 
and the calibration factor a is used to calculate the effective mass m e s of the resonator. However, when 02 ^ there 
is this extra force which is not included in the fitting formula, thus a obtained from the fitting is not equal to oi as 
can be seen in the top plots of Fig. [T7J (01 is chosen to be 10 4 ), and m c ff calculated from a would be incorrect. For 
example, the 0:2/01 = 5 and 10 curves show that o obtained from the fitting is smaller than oi = 10 4 . Since m e g is 
inversely proportional to o, a smaller o would result in a larger m c ff . This is in agreement with our observation that 
the effective mass obtained from electrostatic calibration measurements is larger than expected. The inset of the top 
right plot in Fig. [T7J shows the calibration factor o versus the distance of the closest point used in the fitting with 
experimental data from our electrostatic calibration measurements, showing qualitative agreement with the model. 

The previous analysis indicates that with a carefully chosen combination of p and 02/01, the existence of an extra 
electric force could well explain the problems we experienced in our electrostatic calibration measurements. Another 
very important result from this analysis is that unreliable values of do could be obtained from the fitting if there 
exist forces other than the expected electrostatic force. As shown in the bottom plots of Fig. [17j positive values of 



dp' were obtained, both when the exponent q is left as a free parameter or fixed at 2.5. Based on the way the data 
were constructed, the correct value for do is nm, and the difference between A 



do and its corresponding null value 
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FIG. 16: (Color online) Optimal exponent of the hypothetical scenario data versus the distance of the closest data point used 
in the fitting. The data are constructed from Eq. (|12[) with p = 5, ai = 10 4 , and 0:2/0:1 = 5, 10, 50, 100 for the black circle 
points, red square points, green diamond points and blue triangle points, respectively. 
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FIG. 17: (Color online) Plots of the calibration factor a (top) and contact distance d ' (bottom) of the hypothetical scenario 
data versus the distance of the closest data point used in the fitting. The plots on the left are the results of the best fit when 
the exponent q is a free parameter, and the plots on the right are the results when q is fixed at 2.5. The data are constructed 
from Eq. fJ2]with p = 5, a.\ = 10000, and 02/01 = 5, 10, 50, 100 for the black, red, green and blue curve, respectively. The 
insets on the right are the corresponding results using experimental data from our electrostatic calibration measurements. The 
relationship between the distance and the PZT voltage is d — do — /3(Vp Z T — Vpzt) where /3 is the actuation coefficient of the 
piezoelectric transducer. The same labeling for the points as in the previous figure is used. 



cannot be covered by the fitting uncertainty. This results in a further source of systematic error in the determination 
of the absolute distance, which adds up to other sources, including the one recently discussed in (76J. It should also 
be noted that there seems to be a strong correlation between the two fitting parameters a and rig, which is indeed 
also present in our electrostatic calibration measurements as shown in the insets in the bottom right plot of Fig. IT71 



D. Electrostatic patch effects 



As it is well known, all metallic surfaces in reality are not equipotential surfaces, showing instead voltage variations 
of order 10-100 mV over micrometer distances. These patch potentials are typically due to local changes in the 
work function associated to different crystallographic facets of the metal. Electrostatic patches are known to be an 
important systematic in several precision measurement experiments, including those aiming at detecting the Casimir 
force. Apart from the static component, patch potentials can also fluctuate in time, a dynamical process that has not 
been studied in detail to date. It has recently been shown that by cooling a Au sample the electric-field noise above 



the metal is substantially reduced, a process possibly due to thermal activation barriers in the surface potential [77 1. 

Here we briefly describe the physics of electrostatic patches in the cylinder-plane geometry, and discuss whether they 
could be partly responsible for the anomalous exponent in our electrostatic calibrations. Our considerations follow 
closely the model and notation of (23[. A related effect is the fluctuation-induced interaction between monopolar 
charge disorder within the dielectric slabs [24|. The electrostatic interaction energy between two parallel plates, 
whose surfaces contain stochastic voltage variations V a (x,y) (a denotes the upper or lower plate) is 



— / dtk^—, , 
2 J sinh(fea) 



U PP = ^ I dk ^ un ^ S(k). (14) 



This expression results from assuming zero-average patches, and an isotropic two-point correlation in the transverse 
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plane-wave basis k given by (V^k, V&,k') = S a ^C a ,k5 2 (^ — k'). Here (. . .) denotes stochastic average, k = |k|, and the 
power spectral density S(k) is defined as J °° dkkS(k) = (l/8n) / °° dkk(C\ t k + C 2 .fc). The corresponding electrostatic 
force due to these potential patches is given by F pp = —dU pp /dd. To compute the patch effect on the force in 
the cylinder-plane configuration we make use of the PFA to treat the curvature of the cylinder, which is a good 
approximation in the limit d/a 1. We do not impose any restriction on the typical size of the patches (i.e. we leave 
kd arbitrary). In the limit d/a <C 1 the electrostatic force due to patches in the cylinder-plane configuration is 

p TreoL (d\ 1/2 f°° k^e-^ 

Fcp= vf a UJ Jo dk ^m s{k) - (15) 

Two simple limiting cases can be analyzed. In the large patch limit (kd <C 1) the force is given by 



cp ~ 2V2 d?n V ™-< 



with V 2 ms = J °° dkkS(k). This expression is exactly equivalent to the r.h.s. of Eq. (JTJ) with V 2 replaced by V 2 msl 
as expected for large patches. In the small patch limit (kd >• 1) the force is exponentially suppressed because the 
patches are small and change sign rapidly, resulting in a vanishing net interaction between the plates. It should be 
noticed that the patch force depends on distance as an inverse power law (with exponent 1.5) only in the small-patch 
limit. The smaller the patches, the faster the decay (bigger exponent). More importantly for our purposes is the fact 
that the electrostatic patch force (|16|) is independent of the applied voltage V between the cylinder and the plane. 
Therefore, it cannot explain the anomalous exponent of the electrostatic calibration, that stems from the ^-dependent 
contribution to the force. The patch force ([lo]) is, instead, a background force that could possibly show up in the 
analysis of the electrostatic residuals, that is, the force after the subtraction of the Coulomb-like ^-dependent terms. 

It could be argued that the presence of strong electric fields in the gap between the plates may, in principle, redis- 
tribute the spatial configuration of the patches, and then the question is whether the force between the redistributed 
patches depends on V 2 (note that if this were the case, there might be some hope that the anomalous exponent is 
partially due to patches). Equation (|16p was obtained assuming that the two plates had only stochastic potentials 
fluctuating around zero. If an external fixed potential difference V is applied between the plates, the linearity of 
Laplace's equation implies that the total force will be the sum of the usual V 2 term plus the ^-independent term 
given in (|TB)) . In principle, for sufficiently large external fields in the gap, the power spectrum S(k) could depend on 
the external voltage V, but it is unclear if, and how, such an effect can account for the anomalous exponent in our 
electrostatic calibrations, at least in the ones performed by maintaining the electric-field approximately constant in 
the explored distance range. 



V. MINIMIZING POTENTIAL IN THE PARALLEL PLATES GEOMETRY 



We have also performed in the same experimental conditions measurements in the plane-plane configuration using 
flat coated mirrors facing the resonator. Considering that the sphere-plane geometry has been the subject of former 
work [38l . l39j , this allows us to complete the picture on the relationship between the distance-dependent minimizing 
potential and the specific geometrical configuration. Three mirrors with different coatings (Au, Ag and Al) were used 
to also investigate the possible relationship between Vb and the nature of the substrate. Although the parallelization 
is obtainable with a good level of approximation in one direction only, we can introduce an off-line parallelization 
correction in the data fitting. We analyze separately in the following results on the three substrates. The data of Vq 
from four runs were merged in the same plot for each of the three mirrors as shown in Fig. [TS] where they appear 
versus the curvature coefficient K e \ instead of the distance. This is because the range of K \ achieved is not very large 
especially in the case of aluminum mirrors due to limited parallelization and smaller conductivity from oxide layers 
on the surface of the mirror, and consequently the usual fitting procedure to find the absolute distance would produce 
rather unstable and inaccurate results. Since there is a one to one mapping between K e \ and the distance, K e \ is used 
here as a fair indicator of the distance for all runs. 

In the case of Au mirrors, Vq fluctuates around 100 mV without any noticeable trend. The fluctuation is bigger at 
large distance partly because the frequency shifts tend to be small also for relatively larger bias voltage, this being 
reflected in the error bars in the fitting procedure. As the distance gets smaller, Vo seems to converge to a constant 
value. The minimizing potential Vq shows a similar behavior in all runs, and the values of Vq are also comparable, 
with runs 1, 2, and 4 all around 0.1 V, and run 3 slightly lower at around 0.4 V. In Table III, the average value and 
the standard deviation of Vq for each run is obtained using the six data points with largest K c \ value from each run. 
We see that the standard deviation is reasonably small indicating that Vq is mostly constant. 




FIG. 18: (Color online) Plots of the minimizing potential Vo versus curvature coefficient K e \ from various runs in plane-plane 
geometry with three different coatings of the mirror, Au (left), Ag (middle), and Al (right), all facing the Au-coated resonator. 



Using Ag mirrors instead, a noticeable difference evident from the electrostatic measurement is that we achieve 
values of K e \ smaller than in the case of Au mirrors. This is possibly due to the fact that Ag is easily oxidized once 
exposed in air prior to the insertion into the vacuum chamber. Of course, we cannot rule out the possibility that 
this could also be partly due to slightly different parallclization configurations. In terms of residual potential Vq, the 
Ag mirror shares the same behavior as gold mirror, as Vo fluctuates randomly remaining constant especially at small 
separation gap. However, in this case the Vb value is noticeably higher that that of the Au mirror. 

Finally, with an Al coating, the lowest value of K c \ achieved is much smaller than that of the Au and Ag mirrors. 
This may be due to the easiness to form oxide layers, and it could also be related to the smaller conductivity with 
respect to Al and Au. In this case the residual potential Vb is also constant at small separation gaps as it can be seen 
in Fig. [TS1 and it is manifest that Vb is significantly larger than that for Au and Ag mirrors. A more cjuantitative 
analysis is precluded by the strong dependence of the work function of Al on the exposure time in air [78[. 

A comparison among the various minimizing potentials Vb and physical parameters of the substrate is shown in the 
second part of Table III. The average values of Vb are obtained by averaging three runs from each mirror weighted by 
their variance. From the table it is clear that Vb for mirrors coated with Au, Ag, and Al are significantly different. 
For a rough theoretical comparison, the work function and the Fermi energy as well as their difference are also listed 
in the table. The work function of a metal is closely related to its Fermi energy, but due to the presence of defects 
and impurities on the surface these two quantities do not coincide, with their difference largely due to the surface 
charge distribution and surface dipole distribution. From Table III we can see that Au has the smallest difference 
between its work function and Fermi energy while Al has the largest. Since this difference indicates the magnitude 
of the surface charge and dipole distribution which may be directly related to the residual potential Vb, it could be 
used to explain the different value of Vb measured using mirror with different coatings. In fact the average Vb values 
for Au, Ag, and Al mirrors is consistent with the order of this difference. Regarding the best fitting exponent for the 
scaling of the curvature coefficient with distance, we have obtained the expected one from Coulomb force, although 
it should be noticed that the range of distances was limited to a minimum value of about 3 /xm. 





Run 


Vb (V) 


ov {V) 


<Vo) (V) 


W (eV) 


E F (eV) 


E F - W (eV) 


Gold 


1 


0.113 


0.025 


0.048 


5.1-5.47 


5.53 


0.06-0.43 




2 


0.089 


0.017 












3 


0.041 


0.006 












4 


0.113 


0.008 










Silver 


1 


0.145 


0.017 


0.194 


4.52-4.74 


5.49 


0.75-0.97 




2 


0.182 


0.044 












3 


0.264 


0.044 












4 


0.264 


0.023 










Aluminum 


1 


1.177 


0.033 


1.151 


4.06-4.26 


11.7 


7.44-7.64 




2 


1.169 


0.097 












3 


1.118 


0.036 











TABLE III: Average measured value of Vo for mirrors made of various substrates, weighted average value of Vo, and comparison 
with tabulated data for the work function W and the Fermi energy Ef- 
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VI. ON THE IMPORTANCE OF THE MEASUREMENT AND THE MODELIZATION OF THE 
MINIMIZING POTENTIAL IN CASIMIR FORCE MEASUREMENTS 

Our emphasis on measuring the minimizing potential at all the explored distances for the sphere-plane configuration 
[38l |39| . and for the cylinder- plane and parallel planes configurations described here, is due to the fact that the 
understanding of the minimizing potential dependence on distance and time is crucial for the assessment of the 
accuracy of the Casimir force measurements at small (below ~ 1/irn) distances, and for determining the thermal 
contribution at large (in the 1-5 fxm range) distances. In Table IV we report the formulas for the ideal (perfect 
reflectors, zero temperature) Casimir force and for the Coulomb force in the three different geometries. Following 
[79| , let us define the equivalent Casimir voltage as the external bias voltage that can simulate the ideal Casimir force. 
Due to the different scalings with distance of Casimir and Coulomb forces, this equivalent Casimir voltage must be 
specified at each distance. It is, however, important to point out that the difference between the various geometries 
is just a numerical factor which makes this equivalent voltage for the sphere-plane configuration about half the value 
of the parallel plane case (with the cylinder-plane, as is customary, in between the two extreme cases even from this 
point of view) . The formula for the equivalent Casimir voltage is 



•aw-UJ UJ a- (17) 

where £=360, 192, and 120 for the sphere-plane, cylinder-plane, and parallel plate configurations, respectively. 

In the numerical example presented in the last raw of Table IV, at a distance of 1 /im, which can be considered the 
borderline between the short-distance regime and the long-distance regime in which the thermal contribution starts to 
play a significant role, this equivalent Casimir voltage ranges between 10 and 17 mV depending on the geometry At 
3 /im, where the thermal contribution is expected to contribute as 10-20 % of the total force signal, with an absolute 
value still large enough to be detectable in various apparata, the equivalent Casimir voltage is three times smaller 
(i.e. between 3 and 6 mV). This voltage is of the same order of magnitude of the variation of the minimizing potential 
in a range of few micrometers. To take into account this contribution to properly subtract it from the data, it is 
therefore necessary to model the minimizing potential at the few percent level accuracy. Unfortunately, such stringent 
theoretical characterization of the minimizing potential is not yet available. 

A further degree of uncertainty is also related to the fact that, as pointed out in (39l |. the minimizing potential may 
depend on time, a fact that, for instance, could be attributed to temperature drifts. It has been experimentally shown 
in [80j | using a heated atomic force microscope tip that the contact potential depends on temperature, with a slope 
estimated to be of the order of 4 mV/°C. To perform high-precision tests of the Casimir force, one therefore needs a 
stringent temperature stability of the apparatus during the entire measurement run. In [43| , the observation of fast 
changes in the contact potential have been conjectured as due to the effect of background cosmic rays impinging on 
the apparatus. For a release of about 10~ n C/cm 2 through ionization by cosmic rays at sea level, and considering 
the small values of the capacitances (order of hundreds pF), sudden changes of order 0.1-1 mV could be expected 
(see [8l| for a related discussion). A careful control of the minimizing potential will then require also surrounding 
the apparatus with lead shields to reduce the radiation background or particle detectors to veto the apparatus during 
large ionization events, especially for apparata using microspheres. 





Sphere-Plane 


Cylinder-Plane 


Parallel Planes 


Casimir 


t 3 fc„ R 

360 ' lC d? 




t 2 fjr- S 
240 ~cF 


Coulomb 


-ofV 2 


ire La 1 / 2 T^2 

2V2 d 3 / 2 v 


•so S t/2 


VclM 


/ 2 \ V 2 / j \ V 2 l 
\3(io) \~) d 


[192) {—) d 


/ 2 \ V 2 / X 1 / 2 ! 

{120) \7a~) d 


VSL (1 A™) 


9.85 mV 


13.5 mV 


17.1 mV 



TABLE IV: Summary of relevant formulas for the ideal Casimir force and the Coulomb force in the cases of the sphere-plane, 
cylinder-plane, and parallel plane geometries, with both forces in the first two geometries evaluated using PFA. In the third 
row the equivalent Casimir voltage, i.e. the voltage which needs to be applied in order to simulate the Casimir force at a given 
distance d, is reported. In the last row the concrete value of the equivalent Casimir voltage is reported in the case of a typical 
gap distance of 1 fxm. 
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VII. CONCLUSION 



In this paper we have summarized the main outcomes from our effort to measure the Casimir force in the cylinder- 
plane configuration. The presence of uncontrollable frequency shifts of electric origin at the smallest explored gaps, 
evidenced both from the analysis of the residuals of the electrostatic calibrations and a fast-approach technique, 
prevent us from identifying a Casimir-likc contribution at small distance. At large distances thermal drifts arc large, 
if compared to the expected Casimir force, and a careful control of the dependence of the minimizing potential on 
distance is also required to extract meaningful information about the Casimir force and its thermal corrections. While 
unsuccessful, our search for the Casimir force in this geometry, apart from the development of some data taking and 
analysis techniques applicable elsewhere, has evidenced a number of features which may be of more general interest, 
as we try to summarize in the following. 

First, we have observed anomalous behavior for the best fitting exponent with which the electrostatic coefficient 
is scaling as a function of the cylinder-plane separation. The exponent is significantly smaller than the Coulombian 
one at small distances, while it retains its expected value at the largest explored gaps. In the case of the sphere- 
plane measurements, the exponent was slightly smaller than the expected value but it retained its value in the entire 
explored range of distances, this last being smaller than in the cylinder-plane case due to the smaller electrostatic 
signal available in the sphere-plane configuration [38j . 

Second, we have observed a dependence of the minimizing potential on the cylinder-plane distance, similarly to 
the sphere-plane case. Although its dependence is milder at small distances, it still retains a strong dependence on 
distance at larger gaps, and this requires a careful modeling to subtract its contribution when studying forces in the 
1-5 fim range of interest to discriminate among the various models proposed to incorporate the thermal contribution. 

Finally, we have also explored the case of flat surfaces with a rough parallclization and we have found that in this case 
no anomalous behavior is observed for both the scaling exponent and the minimizing potential, even using different 
substrates for the surfaces. The range of explored distances is definitely limited by the approximate two-dimensional 
parallelism achievable with our setup, and it is therefore unclear if anomalous scaling could be instead observed as in 
the case of the smaller explored gaps in the sphere- plane and cylinder-plane configurations. Both apparata built to 
study the Casimir forces in a parallel plane configuration have not observed dependence of the minimizing potential 
on distance [!, [64| . 

Our findings should be then related to recent outcomes from various experiments confirming the presence of nontriv- 
ial, formerly unidentified systematic effects in the electrostatic calibrations [H, H^, |4l], \^ . Several recent experiments 
are also showing that the observation of Casimir or Casimir-Polder forces is less trivial than previously stated, for 
instance, with regards to the dependence on the optical properties of the substrates [13] and the presence of dielectric 
layers on the substrates (83j . Theoretical arguments have been recently provided for the nontrivial interplay between 
thermal fluctuations and geometry [84[ , thermal, conductivity, and roughness corrections [85|, [8(| and the role of the 
statistical properties of the conducting surfaces [24l ] . Deviations from the pure Coulombian contribution and from the 
hypothesis of a constant minimi zing po tential have also been observed in atomic force microscopy for sharp tips, for 
instance, due to capillary forces [871 188(| . 
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(Argentina). 



[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 (1948). 
[2] M. J. Sparnaay, Physica 24, 751 (1958). 

[3] G. Bressi, G. Carugno, R. Onofrio, and G. Ruoso, Phys. Rev. Lett. 88, 041804 (2002). 

[4] P. H. G. M. van Blokland and J. T. G. Oveerbeek, J. Chem. Soc. Faraday Trans. I 74, 2637 (1978). 

[5] S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 

[6] U. Mohideen and A. Roy, Phys. Rev. Lett. 81, 4549 (1998); B. W. Harris, F. Chen, and U. Mohideen, Phys. Rev. A 62, 
052109 (2000). 

[7] H. B. Chan, V. A. Aksyuk, R. N. Kleiman, D. J. Bishop, and F. Capasso, Phys. Rev. Lett. 87, 211801 (2001); D. Iannuzzi, 

I. Gelfand, M. Lisanti, and F. Capasso, Proc. Nat. Ac. Sci. USA 101, 4019 (2004). 
[8] R.S. Decca, D. Lopez, E. Fischbach, and D. E. Krause, Phys. Rev. Lett. 91, 050402 (2003); R. S. Decca et al, Phys. Rev. 

Lett. 94, 240401 (2005). 




Acknowledgments 



20 



[9] R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepanenko, Ann. Phys. (N.Y.) 
318, 37 (2005). 

T. Ederth, Phys. Rev. A 62, 062104 (2000). 

S. Dimopoulos and G. F. Giudice, Phys. Lett. B 379, 105 (1996). 

Y. Fujii, Nature 234, 5 (1971); Ann. Phys. 69, 494 (1972); Phys. Rev. D 9 874 (1974); Int. J. Mod. Phys. A 6, 3505 (1991). 
E. Fischbach and C. L. Talmadge, The Search for N on- Newtonian Gravity (AlP/Springer-Verlag, New York, 1999). 
V. A. Kuzmin, I. I. Tkachev, and M. E. Shaposhnikov, Pis'ma Zh. Eksp. Teor. Fiz. 36, 49 (1982) [JETP Lett. 36, 59 
(1982)]. 

V. M. Mostepanenko and I. Yu. Sokolov, Phys. Lett. A 125, 405 (1987). 

A. Lambrecht and S. Reynaud, Eur. Phys. J. D 8, 309 (2000). 
R. Onofrio, New J. Phys. 8, 237 (2006). 

B. V. Derjaguin and I. I. Abrikosova, Sov. Phys. JETP 3, 819 (1957). 
J. Blocki, J. Randrup, W.J. Swiatecki, and F. Tsang, Ann. Phys. 105, 427 (1977). 
N. A. Burnham, R. J. Colton, and H. M. Pollock, Phys. Rev. Lett. 69, 144 (1992). 

B. C. Stipe, H. J. Mamin, T. D. Stowe, T. W. Kenny, and D. Rugar, Phys. Rev. Lett. 87, 096801 (2001). 

C. C. Speake and C. Trenkel, Phys. Rev. Lett. 90, 160403 (2003). 

W. J. Kim, A. O. Sushkov, D. A. R. Dalvit, and S. K. Lamoreaux, Phys. Rev. A 81, 022505 (2010). 

A. Naji, D. S. Dean, J. Sarabadani, R. R. Horgan, and R. Podgornik, Phys. Rev. Lett. 104, 060601 (2010). 
J. Mehra, Physica 37, 145 (1967). 

L. S. Brown and G. J. Maclay, Phys. Rev. 184, 1272 (1969). 
M. Bostrom and B. E. Sernelius, Eur. Phys. J. D 8, 309 (2000). 

D. E. Krause, R. S. Decca, D. Lopez, and E. Fischbach, Phys. Rev. Lett. 98, 050403 (2007). 

E. Fischbach, D. E. Krause, V. M. Mostepanenko, and M. Novello, Phys. Rev. D 64, 075010 (2001). 
R. S. Decca, D. Lopez, E. Fischbach, and D. E. Krause, Phys. Rev. Lett. 91, 050402 (2003). 

G. L. Klimchitskaya, R. S. Decca, E. Fischbach, D. E. Krause, D. Lopez, and V. M. Mostepanenko, Int. J. Mod. Phys. A 
20, 2205 (2005). 

R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepanenko, Phys. Rev. D 75, 
077101 (2007) 

R. S. Decca, D. Lopez, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and V. M. Mostepanenko, Eur. Phys. J. C 51, 
963 (2007). 

V. M. Mostepanenko, R. S. Decca, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, and D. Lopez, J. Phys. A 41, 164054 
(2008). 

W. J. Kim, M. Brown-Hayes, D. A. R. Dalvit, J. H. Brownell, and R. Onofrio, Phys. Rev. A 79, 026102 (2009). 
R. S. Decca, E. Fischbach, G. L. Klimchitskaya,D. E. Krause, D. Lopez, and V. M. Mostepanenko, Phys. Rev. D 79, 
124021 (2009). 

D. A. R. Dalvit and R. Onofrio, Phys. Rev. D 80, 064025 (2009). 

W. J. Kim, M. Brown-Hayes, D. A. R. Dalvit, J. H. Brownell, and R. Onofrio, Phys. Rev. A 78, 020101(R) (2008). 
W. J. Kim, M. Brown-Hayes, D. A. R. Dalvit, J. H. Brownell, and R. Onofrio, J. Phys. Conf. Ser. 161, 012004 (2009). 
R. S. Decca, E. Fischbach, G. L. Klimchitskaya, D. E. Krause, D. Lopez, U. Mohideen, and V. M. Mostepanenko, Phys. 
Rev. A 79, 026101 (2009). 

S. de Man, K. Heeck, and D. Iannuzzi, Phys. Rev. A 79, 024102 (2009). 

W. J. Kim, A. O. Sushkov, D. A. R. Dalvit, and S. K. Lamoreaux, Phys. Rev. Lett. 103, 060401 (2009). 
S. E. Pollack, S. Schlamminger, and J. H. Gundlach, Phys. Rev. Lett. 101, 071101 (2008). 
D. K. Schroder, Meas. Sci. Technol. 12, R16 (2001). 

J. B. Camp, T. W. Darling, and R. E. Brown, J. Appl. Phys. 69, 7126 (1991). 
J. B. Camp, T. W. Darling, and R. E. Brown, J. Appl. Phys. 71, 783 (1992). 
J. B. Camp and R. B. Schwarz, Appl. Phys. Lett. 63, 455 (1993). 

M. Bordag, B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. Lett. 85, 503 (2000). 
C. Genet, A. Lambrecht, and S. Reynaud, Phys. Rev. A 62, 012110 (2000). 
C. Genet, A. Lambrecht, and S. Reynaud, Int. J. Mod. Phys. A 17, 761 (2002). 

B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. A 65, 062109 (2002). 
V. B. Svetovoy and M. V. Lokhanin, Phys. Rev. A 67, 022113 (2003). 
J. S. H0ye, I. Brevik, J. B. Aarseth, and K. A. Milton, Phys. Rev. E 67, 056116 (2003). 
B. Geyer, G. L. Klimchitskaya, and V. M. Mostepanenko, Phys. Rev. A 67, 062102 (2003). 
R. Esquivel, C. Villarreal, and W. L. Mochan, Phys. Rev. A 68, 052103 (2003). 
J. R. Torgerson and S. K. Lamoreaux, Phys. Rev. E 70, 047102 (2004). 
R. Esquivel and V. B. Svetovoy, Phys. Rev. A 69, 062102 (2004). 

I. Brevik, J. B. Aarseth, J. S. H0ye, and K. A. Milton, Phys. Rev. E 71, 056101 (2005). 
V. S. Bentsen, R. Herikstad, S. Skriudalen, I. Brevik, and J. S. H0ye, J. Phys. A 38, 9575 (2005). 
S. K. Lamoreaux and W. T. Buttler, Phys. Rev. E 71, 036109 (2005). 

A. Lambrecht, V. V. Nesvizhevsky, R. Onofrio, and S. Reynaud, Class. Quantum Grav. 22, 5397 (2005). 
P. Antonini, G. Bressi, G. Carugno, G. Galeazzi, G. Messineo, and G. Ruoso, New. J. Phys. 8, 239 (2006). 
G. Bressi, G. Carugno, A. Galvani, R. Onofrio, G. Ruoso, and F. Veronese, Class. Quantum Grav. 18, 3943 (2001). 
P. Antonini, G. Bimonte, G. Bressi, G. Carugno, G. Galeazzi, G. Messineo, and G. Ruoso, J. Phys. Conf. Ser. 161, 012006 



21 



(2009). 

D. A. R. Dalvit, F. C. Lombardo, F. D. Mazzitelli, and R. Onofrio, Europhys. Lett. 67, 517 (2004). 

M. Brown-Hayes, D. A. R. Dalvit, F. D. Mazzitelli, W. J. Kim, and R. Onofrio, Phys. Rev. A 72, 052102 (2005). 

M. Brown-Hayes, J. H. Brownell, D. A. R. Dalvit, W.-J. Kim, A. Lambrecht, F. C. Lombardo, F. D. Mazzitelli, S. M. 

Middleman, V. V. Nesvizhevsky, and R. Onofrio, J. Phys. A 39, 6195 (2006). 

T. Emig, R. L. Jaffe, M. Kardar, and A. Scardicchio, Phys. Rev. Lett. 96, 080403 (2006). 

M. Bordag, Phys. Rev. D 73, 125018 (2006). 

H. Gies, K. Langfeld, and L. Moyaerts, J. High Energy Phys. 06 (2003) 018. 

H. Gies and K. Klingmuller, Phys. Rev. Lett. 96, 220401 (2006); Phys. Rev. D 74, 045002 (2006). 

W.J. Kim, M. Brown-Hayes, D. A. R. Dalvit, and R. Onofrio, in Proceedings of the XLIInd Rencontres de Moriond, 
Gravitational waves and experimental gravity, J. Dumarchez and J. Tran Thanh Van editors (The Gioi Publishers, 2007), 
pp. 291-298. 

W. R. Smythe, Static and Dynamic Electricity (McGraw-Hill, New York, 1968). 
H. Gies and K. Klingmuller, Phys. Rev. Lett. 97, 220405 (2006). 
Q. Wei and R. Onofrio, Phys. Lett. A 374, 2230 (2010). 

P. J. van Zwol, V. B. Svetovoy, and G. Palasantzas, Phys. Rev. B 80, 235401 (2009). 

J. Labaziewicz, Y. Ge, D. R. Leibrandt, S. X. Wang, R. Shewmon, and I. L. Chuang, Phys. Rev. Lett. 101, 180602 (2008). 
M. Uda, Y. Nakagawa, T. Yamamoto, M. Kawasaki, A. Nakamura, T. Saito, and K. Hirose, J. El. Spectr. 88, 767 (1998). 
R. Onofrio and G. Carugno, Phys. Lett. A 198, 365 (1995). 

J. Remmert, Y. Wu, J. Lee, M. A. Shannon, and W. P. King, Appl. Phys. Lett. 91, 143111 (2007). 

V. P. Mitrofanov, L. G. Prokhorov, and K. V. Tokmakov, Phys. Lett. A 300, 370 (2002). 

V. B. Svetovoy, P. J. van Zwol, G. Palasantzas, and J. Th. M. De Hosson, Phys. Rev. B 77, 035439 (2008). 

G. Palasantzas, V. B. Svetovoy, and P. J. van Zwol, Phys. Rev. B 79, 235434 (2009). 

A. Weber and H. Gies, Phys. Rev. D 80, 065033 (2009). 

P. A. Maia Neto, A. Lambrecht, and S. Reynaud, Europhys. Lett. 69, 924 (2005). 

A. Canaguier-Durand, P. A. Maia Neto, I. Cavero-Pelaez, A. Lambrecht, and S. Reynaud, Phys. Rev. Lett. 102, 230404 
(2009). 

J. Jang, G. C. Schatz, and M. A. Ratner, J. Chem. Phys. 120, 1157 (2004). 

G. Huber, H. Mantz, R. Spolenak, K. Mecke, K. Jacobs, S.N. Gorb, and E. Arzt, Proc. Natl. Acad. Sci. USA 102, 16293 
(2005). 



